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Abstract 

In this work, nonreflecting boundary conditions in generalized three-dimensional curvilin¬ 
ear coordinates are derived, relying on the original analysis that was done in Cartesian two- 
dimensional coordinates by Giles jl] . A thorough Fourier analysis of the linearized Euler equa¬ 
tion is performed to determine the eigenvalues and the eigenvectors that are then used to derive 
the appropriate inflow and outflow boundary conditions. The analysis lacks rigorous proof of 
the well-posedness in the general case, which is open to investigation (a weak assumption is 
introduced here to complete the boundary conditions). The boundary conditions derived here 
are not tested on specific applications. 


1 3D Linearized Euler Equations in Cartesian Coordinates 

The non-linear, non-conservative (primitive) form of the Euler equations in Cartesian coordinates 
can be compactly written: 


Qt + + BQy + — 0 
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where Q = { p u v 
matrices. 


w p is the vector of the primitive variables and A, B and C are 5x5 
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Here, p is the density, p is the pressure, and m, v and w are the velocity components. The 
linearization is done by decomposing the primitive vector Q into a mean part and a perturbation 
part: 


Q — Q + Q 


( 3 ) 


where 


Q=(p u V w p and Q.' = { p' 


u V w 


> / 
j p ) 


( 4 ) 


and because the perturbations are assumed to be much smaller than the mean flow, only the terms 
of order O(Q') are kept, while the other terms of order 0(Q^^) are neglected. Thus the linearized 
Euler equations, after the variables are non-dimensionalized using the mean density and the speed 
of sound, are: 


q* + aq; + bq; + cq; = o (s) 

where the matrices A, B and C are defined as 
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2 3D Linearized Euler Equations in Curvilinear Coordi¬ 
nates 

Consider the transformation between the curvilinear coordinates and the Cartesian coordinates: 


T = rit) 

V = v{x,y,z,t) (7) 

C = C{x,y,z,t) 

Assume for now that the grid is not moving: the time dependence in the above relations is 
dropped out. Using the chain-rule formulation applied to Eq. ([5]), the linearized Euler equations are 
written in curvilinear coordinates system as: 
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( 8 ) 


q; + [e.A + ^B + e.^Q' 

+ [t]^A + T]yB + r]^c] Qy 

+ [CxA + (^yB + = 0 

Now denote 


A — ^a;A + ^j,B + 

B = r]^A + T]yB + r]^C 
C = C.A + CyB + CzC 

As a result, Eq. ([8]) can be written as: 


Qj + AQ^ + BQj^ + CQ^ — 0 
The matrices A, B and C are: 
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The contravariant mean velocities are given by: 


U = ixU + iyV + 

V = rjjl + rjyV + rizW (14) 

W = C^U + CyV + CzW 

3 Fourier Analysis of the 3D Linearized Euler Equations in 
Curvilinear Coordinates 

A three-dimensional wave-like solution in the form: 


Q' _ Q^iiM+l-n+mC-ujt) 


(15) 


is considered and introduced into Eq. flTOl) . This will generate a matrix equation for the 3D wave 
amplitudes as unknowns. The matrix has the form: 
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(18) 








Eq. (fT6|) or ffT7|) is a homogeneous set of equations; in order to admit non-trivial solution the 
determinant of the matrix must satisfy: 


det{—ujl + kA + /B + mC) = 0 


(19) 


or 


(3^ (/5^ — al — al — al) = 0 


( 20 ) 


3.1 Eigenvalues 

If k is the unknown of Eq. fl20|) the hrst three roots are identical: 

oj — VI — Wm 

h-\ 2 3 — -^- 

u 

The other two roots are determined by solving the equation: 

{Uk + Vl + Wm — — {^xk + rjxl + 

- {iyk + r]yl + - {izk + rjV + = 0 

Denote: 


lel 

Ihl 

ICI 

\id\ 

m 

IhCI 

h 

T 

The discriminant is: 
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( 21 ) 
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So the fourth and the hfth roots are: 


or 


^4,5 — 


—2 (S + /rt/) ± \/A 

2(|{|2-r) 
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^ 4,5 


(S + /iI7) {-l±S) 


(24) 


s 



(kP - U^) (T - 1,^) 

(2 + l^uf 


(25) 


If u is complex with Im{u) > 0 then the right propagating waves are those for which Im{k) > 0. 
This is because the amplitude of the waves is proportional to exp[Im{u){t — i/c)] where c is the 
aparent velocity of propagation of the amplitude. If u: and k are real then the standard result is 
the analysis of dispersive wave propagation using the group velocity. For real u the incoming waves 
are those which either have Im{k) > 0, or have real k and (Ufe > 0 To determine the direction of 
propagation of the waves, let’s solve Eq. (120|) for uj too. The roots are easily found: 


1 ^ 1 , 2,3 = Uk + Vl + Wm 


(26) 


= Uk + Vl + Wm ± Jai + + o| (27) 

The analysis of the hrst three waves is straigthforward. If Im{uj) > 0 then it is clear that 
Im{k) > 0. If Im{u) = 0 the group velocities for the hrst three roots are all equal to U which means 
that, for U > 0 the corresponding waves are incoming waves at a; = 0 and outgoing waves at a; = 1. 
For the fourth and the hfth roots the group velocities are: 


Cg4,5 


Jj ^ ‘^0‘l^x + ‘^CX2^y + 2 q : 3(^2 

2 T 0.2 T ttg 


(28) 


If the I and m components of the wavenumber are small (hypothesis used to derive the approx¬ 
imate boundary conditions) the second term of the right-hand-side of Eq. (|28|) tends to |^|. The 
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fourth wave is incoming wave at x = 0 and outgoing wave at x = 1. For 0 < t/ < |^| (subsonic 
case) the hfth wave is an outgoing wave at x = 0 and an incoming wave at x = 1. For supersonic 
case {U > 1^1) the hfth wave behaves as the fourth wave (incoming wave at x = 0 and outgoing 
wave at X = 1). It can be proved that if u and/or S are complex then one of the two roots for 
k has imaginary part, while the other has a negative imaginary part: is dehned to be the root 

with positive imaginary part (right-running wave), and k^ is a root with negative imaginary part 
(left-running wave). 

Remark: U is the contravariant velocity: the relation 0 < < |^| is translated into Cartesian 

coordinates as 0 < u < 1, taking into account that = 1, = 0 and = 0 (in Cartesian 

coordinates). 

Next the right and the left eigenvectors corresponding to matrix equations flTbl) and ffT7)l . respec¬ 
tively, will be determined. 


3.2 Eigenvectors 

3.2.1 Root 1: entropy wave 


uj — VI — Wm — — — 

ki = -=-, u: = Uki + Vl + Wm 


(29) 


The right eigenvectors are determined by solving the matrix equation (ITHll . First, denote 
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The general solution (depending on 3 independent constants Ci, C 2 and C^) to this set of equa¬ 
tions is given by: 
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The choice of these constants will give the hrst three right or left eigenvectors. The hrst right 
eigenvector (keeping the notation of Giles [I]) is chosen by setting C\ = —1/|^| and (^2 = Gs = 0. 
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The left eigenvectors are determined by solving the matrix eqnation flTTI) which becomes: 


(33) 


( h h h U h ) 
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(34) 


The general solntion (depending on 3 independent constants Di, D 2 and D^) to this set of 
eqnations is given by: 


( ^"1 h h h h ) — { -^1 0:2-02 + O3-O3 — 0:1^2 ~ OiiD^ — Di ^ 
The hrst left eigenvector is chosen by setting Di = —1.^|, O 2 = 0 and = 0: 

uf = ( -1^1 0 0 0 1^1 ) 

The vector vf' is obtained: 


lim (k;) 

\i^0,X2^0 
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(35) 


(36) 










(37) 


A = ( -151 0 0 0 151 ) 

This choice of the eigenvector corresponds to the entropy wave traveling downstream at a speed 
U. This can be verihed by noting that the only non-zero term in the right eigenvector is the density, 
so that the wave has varying entropy, no vorticity and constant pressnre. The left eigenvector 
measnres entropy in the sense that u^U is eqnal to the linearized entropy, ^x{p' — p')- The right 
eigenvector uf is chosen snch that the left eigenvector is normal to it. 

3.2.2 Root 2: first vorticity wave 
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^ = Uk 2 + VI + Wm 


(38) 


The second right eigenvector corresponding to the triple-root is chosen by setting Ci = Ca = 0 
and 6*2 = — ^ in Eq. fl3^ . Taking into acconnt that 0.2 = ^yk 2 ^riyl + C,ym and Oi = 
and replacing ^2 using Eq. fl38|l . the right eigenvector is obtained as 
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(39) 


and still keeping a mnltiplicative constant, C 2 , for the reqnirement that the left eigenvector be nor¬ 
mal to the right eigenvector. This choice is not nniqne. Any linear combination of the eigenvectors is 
itself an eigenvector, and the only constraint is the reqnired orthogonality condition. The motivation 
of this choice is a knowledge of the distict behavior of the entropy and vorticity variables in flnid 
dynamics; this choice will also simplihes the algebra at later stages of this analysis. 

The second left eigenvector corresponding to the triple-root is chosen by setting Di = D 3 = 0 
and D 2 = —^ in Eq. (I35|) . Taking into acconnt that 02 = ^yk 2 + riyl + Cym and oi = ^xk 2 + Pxi + Cx’m, 
and replacing ^2 using Eq. (I38|) . the left eigenvector is obtained as 


U 2 — ( 0 /2,2 h,3 00 ) (40) 

where 

12,2 = - [4 (1 - FAl - WA 2 ) + iVyXl + CyX2) U] 

4,3 = (1 — EAi — WX 2 ) + {pxXi + (xX 2 ) (41) 
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Now the constant must be determined such that u^] = 1. Only the right 

eigenvector is modihed, and its new form is: 


where 


u? = 




/ 0 \ 

- [4 (1 - FAi - TFA 2 ) + (r/,Ai + C.As) U] 

— y^i — w A 2 ) + + CxA2) F] 
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The vector is obtained: 


^2 = xlii + 


lim ikX) = = 
Ai^-0,A2^’0 U 


(42) 


(43) 


^2 — =^2 A — ( 0 /2,2 4,3 0 4,5 ) 


(44) 


where 


4,5 — iy + Ca;-^2) “ (hy^l + Cy^2) 

This root correponds to a vorticity wave rotating around and traveling downstream at a 

speed U. 

3.2.3 Root 3: second vorticity wave 
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u — VI — Wm 
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a; = f//c 3 + 1 // + Wm 


(45) 


The _third right eigenvector corresponding to the triple-root is chosen by setting Oi = 0*2 = 0 and 
Os = — ^ in Eq. (132|1 . Taking into account that 03 = + 'Hzl + and oi = F ^3 + Vxl + Cx'm, 

and replacing ks using Eq. (05]), the right eigenvector is obtained as 
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[e. (1 - FAi - TFA2) + (r/.Ai + CxAa) U] 

0 
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(46) 


and still keeping a multiplicative constant, C^, for the requirement that the left eigenvector be 
normal to the right eigenvector. 

The third left eigenvector corresponding to the triple-root is chosen by setting Di = D 2 = 0 and 
D 3 = —— in Eq. (I5S]) . Taking into account that 03 = -f rjJ + Qm and ai = + rjJ + 

and replacing using Eq. fl45|l . the left eigenvector is obtained as 


.3 — ( 0 [3^2 0 /3^4 0 ) 


(47) 


u: 


where 


h,2 — ~ [F (1 ~ — hEA2) -l- -I- CzA2)F] 

4,4 = [F (1 - FAi - IFA2) + (r^.Ai + aA 2 ) U] 


(48) 


The constant must be determined such that [ug u^] = 1. Only the right eigenvector 

is modihed, and its new form is: 
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where 


>i'3 = viJTi 


(50) 


The vector Vg is obtained: 


lim 

Ap —^0, A2 —^0 
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(51) 


where 



A — ( 0 4^2 


0 4,4 



4,5 ~l~ Cx^2^ (hz-^l ~l~ C 2 '^ 2 ) 

This root correponds to a vorticity wave rotating around r^-axis, and traveling downstream at a 
speed U. 

3.2.4 Root 4: first acoustic wave 

Using the notations in Eqs. flTSl) and fl23|) and also the Eq. ([25]); 

A ;4 = ^ ^ + VI + Wm + Ja\ + 0 ^ + ai (53) 

(l«P-Cl) 

The matrix equation used to determined the right eigenvectors for this eigenvalue is: 

/3 oi 02 03 0 

0/5 0 0 oi 

0 0 /5 0 0(2 

0 0 0 /5 Os 

0 CKi 0.2 03 /5 

where 


( fi\ 

h 
h 

f'A 

\h J 


= 0 


(54) 


/5 = f//c 4 + Vl + Wm — u = —^^a\ + al + al 
The general solution (depending on a constant C 4 ) is: 
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The first right eigenvector is obtained by setting 6*4 = -: 


/ 1-Ukl- l^Ai - ITA 2 \ 
^xkl + + CxA2 

^yk\ + + Cy^2 

^zkl + T]z\l + Cz-^2 

V 1 - FA;* - FAi - TFA2 / 

where is an arbitrary multiplicative constant and 


kl = 


{E* + ij*U) i-l + S*) 


with the new dehnitions in terms of Ai and A 2 : 


= i-VX^-WX2 
E* = Ai|er/| + A2|eCI 

T* = Xl\v\^ + Xl\C\^ + 2X,X2 \vC\ 


S* = 
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* + iiHff' 


The left eigenvector is obtain by solving the matrix equation: 


( h h I 3 U h ) 


(3 Oil 0^2 OiS 0 
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which result in the general solution (depending on a constant ^ 4 ) 


= 0 


( 12 is ^4 ^5 ) — (0 aiD^ q;2-D4 asD/i (3Di ) 

The left eigenvector is obtained by setting D/i = 


(57) 


(58) 


(59) 
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U4 = C'f ( 0 /\2 llz 1*4,4 h,5 ) (60) 

where is an arbitrary multiplicative constant and 

^4,2 ~ {ixkX + + Cx-^ 2 ) 

(^4,3 = {iy^4 + + C?/-^2) 

i*4,4 — {izk*4 + ^7zAi + GA 2 ) 

= 1 _ I7A;* - FAi - TTAa 

The constants and Ff must be determined such that [U 4 uf] 
of the eigenvectors u;f and U 4 are: 

l-Ukl- FAi - iFAa A 

+ hxAl + CxA2 

^yk-X + + C 1 /A 2 (62) 

izkX + Tj^Xl + CzA2 
1 - FA;* - FAi - fFA2 / 

and 



(61) 

= 1. The new forms 



(63) 

(64) 

(65) 


14,2 = ^x + \i{UVx-VQ+\ 2 {UC.-W^ 

h,3 = ^y + \l{UVy-ny)+X2{UCy-WQ 

1,4 = F + Ai(Fr/,-FF)+A2(FF-WF) (66) 

/4,5 = F+A;:(|ep-F')+Ai(0^,-FF)+A2(0^c-W) (67) 

This wave corresponds to an isentropic, irrotational acoustic wave, traveling downstream. 
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3.2.5 Root 5: second acoustic wave 

Using the notations in Eqs. ffTSl) and fl23ll and also the Eq. fl25ll : 


^ a; = Ufcs + U/+ lUm - 

(leP-t/) ^ 

The matrix equation for this eigenvalue is: 
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where 

/3 = Uk^ + VI + Wm — 00 = a\ + 0^ + Oi\ 
The general solution (depending on a constant C^) is: 
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The hrst right eigenvector is obtained by setting Us = — 
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~ (Cx^5 + hxAl + Ca:A2) 
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where C§‘ is an arbitrary multiplicative constant and 


(71) 
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,. (S* + 1 ,‘U) (-1 - S-) 

with the new dehnitions given by Eq. fl59p 

Following the same routine, the left eigenvector is obtain as: 




n 1 * 1 * 1 * 1 * 

''5,2 '5,3 '5,4 '5,5 


where where is an arbitrary constant and 


^5,2 — ~ + TjxXl + CxX2) 

^5,3 = - (Cj/^5 + Vy^l + Cy^2) 

^ 5 , 4 : — ~ i^zkl + T]zXl + CzX 2 ) 

= Tjk* + VX4+WX2-l 

The constants and must be determined such that [ug u^] 
of the eigenvectors and Ug are: 
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Taking into account that: 


the vector Vg is obtained: 


lim (fcg) 
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(72) 
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(74) 

1. The new forms 
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(78) 
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u-\^\ ' 


;A = (0 i. 


5,2 


h,3 h,4 



where 

4.2 = -^.-Xi{UVx-V^,)-X2{UC,-W^,) 

4.3 = -^y-Xi{UVy-V^y)-\2{UCy-WQ 

4.4 = -^z-Xi{UVz-VQ-X2 {UC-WQ (79) 

4.5 = -F-A;*(|ep-I7')-Ai(0^^-W)-A2(0^C-'^) 

This wave corresponds to an isentropic, irrotational acoustic wave, traveling upstream. 


4 First-Order Unsteady Boundary Conditions for 3D Curvi¬ 
linear Euler Equations 

The one-dimensional, non-reflecting boundary conditions are obtained by ignoring all variations in 
the y- and z- directions and setting Ai = 0 and A 2 = 0. In these conditions: 


= 0, T* = 0, /i* = 1, S* = = 
^ U 


(80) 


|ai=0,A2=0 


(81) 


^nlAi=0,A2=0 ^nlAi=0,A2=0 


The right eigenvectors w^|„=i 2 , 3 , 4,5 are: 
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Ti 


/ 0 \ 
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V 0 / 


17 



(83) 



The left eigenvectors w^|„=i 2 , 3 , 4,5 are: 

= ( - 1^1 0 0 0 1^1 ) 

W2 = (0 - ^2/ Cx 0 0 ) 

wf = ( 0 -G 0 0 ) (84) 

W 4 = ( 0 1^1 ) 

W5 = (0 -ix -iy -iz 1^1) 

The transformation to and from 1-D characteristics variables is given by the next two matrix 
equations: 



where ci C 2 C 3 C 4 and C 5 are the amplitudes of the hve characteristics waves. Let’s check the consis¬ 
tency with the Cartesian coordinates. If 

e. = 1, 4 = 0, 4 = 0 

Vx = 0, = 1, ?7^ = 0 (87) 

4 = 0, 4 = 0, 4 = 1 
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then the Eqs. flH^ and become: 
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(88) 


(89) 


being the same with the results from the thesis of Shivaji [5]. 

At the inflow boundary, the correct unsteady, non-reflecting boundary conditions for a subsonic 
flow are: 


/ Cl \ 


C2 

C3 


= 0 


V 


C4 


/ 


(90) 


while at the outflow boundary: 


C5 — 0 


(91) 


5 Approximate, Quasi-3D, Unsteady Boundary Conditions 
for Curvilinear Euler Equations 

The left eigenvectors are all functions of Ai and A 2 which are assumed to be small (quasi-3D). It 
is prefered to have these functions in polynomial form; the straightforward way to transform them 
into polynomial functions is by using Taylor series expansions. The expansions are written as: 
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(92) 


(Al, As) = v^(0, 0) + Al ((v^)ai);,^^;,^^o + ^2 ((Vn)A 2 );,^^;,^^o + HOT, 

n = 1,2, 3,4, 5 

where HOT stands for higher-order-terms. Equations (l9^ are exact. 

The hrst order approximation which was discussed in the previous section is obtained by keeping 
only the hrst term in the Taylor series expansion. The second order approximation is obtained by 
taking into account the hrst three terms and neglecting the HOT. Higher order approximation is 
possible by considering more terms in the series, but the algebra becomes more complex. 

The second order approximation is: 


(Al, As) ^ (Al, As) = v^(0, 0) + Ai i{^n)xi) x,=X ,=0 + ^2 ii^n)x2) x,=x ,=0 


Ai=A2=0 

m f k. 


= u^(0,0) + -^ (—(u^)aiA) + T: I 

^ \ ^ / Ai=A2=o ia; \ w / Ai=A2=o 


This produces the boundary conditions: 


(93) 


wu;|a,.a,=o + i I —(uJ')a,A 


UJ 


Ai=A2=0 


+ m I ^(u^)a2A 


Q = 0 


Replacing ca, I and m hy —i-§^ and respectively, Eq. 


Ai=A2=0 

becomes: 


(94) 



kn 

UJ 


>n)AiA 


Qr: 


k-n 

UJ 


u; 


()a2aJ 


= 0 


Ai=A2=0 


(95) 


The next step is to take the derivatives (u^)ai and (u^)a 2 , and evaluate u^, ^^(u^)AiAj and 

^^(u^)A 2 Aj in the hypothesis that Ai = 0 and As = 0. This will give hve equations representing 
the approximate, quasi-3D boundary conditions: the hrst four equations are solved at the inhow, 
and the hfth at the outhow. 

At the inhow, the boundary conditions are: 


■-lel 

0 

0 

0 


0 0 0 
-4 0 
-e. 0 


lel 

0 

0 

lel 




^11 

5'12 

5'13 

1714 



921 

922 

923 

1724 








9'ii 

932 

933 

5'34 



. 9ii 

9i2 

943 

1744 


hii 

hi2 

hl3 

hl4 

^15 


hsi 

hs2 

to 

CO 

hs4 

hs5 


hsi 

^32 

CO 

CO 

^34 

^35 


/l4i 

^42 

^■43 

/i44 

^45 


1715 



5'35 
1745 . 


Qc = o 


(96) 
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and at the outflow, they are: 


where 


;o -e. -4 


I'd ] Qt + [ S'si 

+ [ hsi /l52 


fi'52 9bi 955 ] Q,, 

^53 ^54 /l55 ] = 0 


and 


9u 

912 

913 

fl'M 

915 

921 

922 
923 
92A 
925 

931 

932 

933 
fl'34 
935 
9ai 
9A2 
9A3 
9aa 

9A5 

951 

952 

953 
95A 

955 


0 

0 

0 

0 

0 

0 

UVy - V^_y 

-ut]^ + 

0 

^yVx 

0 

0 

-Ut]^ + V^a; 

^xVz ^zVx 

0 

-Ut]^ + V^a; 

-Urjy + Viy 

-Ug. + V^z 

0 

UVx-Vix 
UVy - ny 
Uv. - 


(97) 


(98) 
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hii = 0 

hi 2 = 0 
hi3 = 0 

/ll4 = 0 

/ll5 = 0 

/l21 = 0 

h22 = UCy-W^ 
h23 = -UC. + W^, 

^24 = 0 

^25 ^xCy ^yCx 

^31 = 0 

h 32 = UC .- Wi , (99) 

^33 = 0 

h3A = -UCx + Wi^ 

^35 ^xCz ^zCx 

/z .41 = 0 

K 2 = -uCx + w^:, 

h,3 = -UCy+My 

/l44 = -UCz+W^, 

^51 = 0 

h2 = UCx-W^, 
h3 = UCy-W^y 
hi = Uh-Wi, 

h, = -fJ^+W|5| 

It can be verified that the above equations reduce to the Cartesian 3D Giles boundary conditions. 
Thus, according to Eq. (|87)l . the Eqs. (|^ and (|7n) become: 


■ -1 

0 

0 

0 

1 ■ 


■ 0 

0 

0 

0 

0 ■ 

0 

0 

1 

0 

0 

Qt + 

0 

u 

V 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

V 

0 

0 

1 

0 

0 

1 _ 


_ 0 

V 

—u 

0 

V _ 
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0 0 0 0 0 

0 0 w 0 0 

0 M 0 uJ 1 

0 wJ 0 —u w 


and 


Q. = o 


[0 -1 0 0 1 ] Qi + [ 0 -V u 0 n ] 

+ [ 0 —w 0 u w ] = 0 


being consistent with the results from the thesis of Shivaji 


( 100 ) 


( 101 ) 


6 Analysis of Well-Posedness 

The one-dimensional approximation of the boundary conditions is always well-posed (demonstrated 
by using the energy method). The second order approximation need to be analyzed and corrected 
in case of ill-posedness. This follows next. 

6.1 Inflow Boundary Conditions from the Second-Order Approximation 

The aim is to verify that there is no incoming mode which exactly satisfies the boundary conditions. 
The analysis is done in a frame of reference that is moving with speed V in 77 -direction and W 
in ^-direction. The well-posedness in this frame of reference is both necessary and sufficient for 
well-posedness in the original frame of reference. This simplifies the algebra by setting V = W = 0. 
If there are N' incoming waves, then the generalized incoming mode may be written as 




,n=l 




with Im{u) > 0. The wavenumbers (divided by frequency) are given by: 


where 


k* — k* 


1 

V 


(S* + U) (-1 + S*) 

(kp-p') 


( 102 ) 


(103) 


(104) 
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(105) 


= 


1 - 




- 1 ) 


\ i^*+W 


and the definitions given by Eq. fl5^ . Following the theory given in Giles work [], 
matrix C[ 4 x 4 ] corresponding to the infiow boundary has the elements: 


Cnj = 


The vectors 


= u„|ai=a2=o +Ai (fc;(u„)AiA) + Aa (fc;(u„)A2A . 

Ai=A2=0 V / Ai=A2=0 

n = 1, 2, 3,4 


are needed first with V = W = 0. They are: 




— ( h,i h,2 h,3 ) 


where 


h,i — 

h,2 = 0 
^1,3 = 0 
h,4 = 0 

^1,5 

V 2 = ( h,l h,2 h,3 ^2,4 h,5 ) 

where 

^ 2,1 = 0 

h,2 = —^y — ^iUr]y — X2UCy 

h,3 = ^xXiUrjx + X 2 UCX 

= 0 

^ 2,5 = Ai {^xVy ^yVx) A 2 {^xCy ^yCx) 


the critical 

(106) 

(107) 

(108) 

(109) 

( 110 ) 

( 111 ) 
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(112) 


~ { h,i h,2 h,3 h,A h,5 ) 

where 

^3,1 = 0 

h,2 = —^z — ^iUr]z — \2U(z 

h,3 = 0 

^ 3,4 = Cx + ^iUt]x + \2U Cx 

^3,5 -^1 i^xVz ^zVx) -^2 (^xCz ^zCx) 




( ^4,1 h,2 U,3 U,A U,5 ) 

where 

^ 4,1 

= 0 


U,2 

= ^x + AiUr]x + \ 2 UCx 


h,3 

= ^y AlUf]y + X 2 U(^y 


h,A 

= ^z + AiUr]z + A2 U(z 


h,5 



Thus, the elements of the critical matrix C are: 


Cll 

Cl2 

Cl3 

Ci4 

C21 

C 22 

C 23 

C 24 

C 3 I 


1 

0 

0 

0 

0 




0 

0 


(^y -h AlUfJy -h A2I7 Cy) + (^x + + A2I7Cx) 

[fe + + A2U(y'j + AiUrjz + A2l7Cz)] 


(113) 


(114) 


(115) 


(116) 
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C32 

C33 

C34 

C41 

C42 

C43 

C44 


1 


1 

0 

0 

0 

0 


[(■C?; + ^iUi]y + \ 2 UCy) (.^2 + XiUrjz + A2^7Cz)] 

[ {^z + XiU rjz + \2U (^2) + (Ca; + ^1^7 Tjx + A2^7 Cx) 


{u + \^\) 
2|ep 


kiej + Xiuel + Ximi + Aiier^l + A2|eci + 2AiA2l7|r/C| ' 

+KU (a,|«,,| + A 2 |«I - 65 - AiI7l|l - A2Z71|1) 

+04 + AiC7M + AjC/M 


The requirement that there is no non-trivial incoming mode satisfying the boundary conditions 
is equivalent to the statement that the determinant of the above matrix is non-zero for all real I and 
m and complex u with Im{u) > 0. 

In other words, the procedure would be this: 

-assume that the determinant det{C) = 0; 

-if after manipulations we get FALSE, then the problem is well-posed 
-if we get a combination of u, I and m that make det{C) = 0, then the problem is ill-posed 
Assuming that everything else (grid metrics, mean flow) is constant, the determinant det{C) is 
a function of Ai and A 2 . So, det{C) = 0 will produce a relation between Ai and A 2 (unless the result 
is FALSE). 

To simplify the analysis a little bit, let’s assume that the grid is orthogonal at the boundary. 
Under this hypothesis: 


where 


l^hl = leci = \vC\ = 0 

(117) 

[USs (04 - IT) + Z 7 r + 84] 

(118) 

K= A-2(-i + y) 

(119) 

ej-u 

^• = =\/(c^'-e|)r + e| 

(120) 
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(121) 


r = Ajej + A=e= 

One of the result of det{C) = 0 is C 44 = 0. Replacing kl will give the equation: 




0 ;= 


U 


Replacing S*, and solving for F, gives: 


r = A?0; + A202 = 


The conclusion is: 


if u = +iU 


satisfying the condition that Im{uj) > 0 , then: 


I12Q2^^2Q2 


e| 


( 122 ) 


(123) 


(124) 


e? 


e? 


A^0^ + A 20 a = —=1 and S* = ^ 


U 


U 


(125) 


with the correct branch of the square root being taken to ensure that Im{ki) > 0 . 

To evaluate the critical matrix under the conditions given by Eqs. 01171) and 01251) . denote first: 


71 = ^ 2 ; + ^lU rjx + \ 2 U C,x 

12 = iy ^lUrjy + X 2 U(y (126) 

73 = 'Cz + Vz + ^ 2 U Cz 


The critical matrix is: 


C = 


1 0 

0 ^( 7 | + 7 ?) 
0 ( 7273 ) 

0 0 


0 0 

^ (7273) 0 

^(7l + 7?) 0 

0 0 


(127) 


Next, the second row is multiplied by 72 , and the third row by 73 , then add the second row to 
the third row. The result is: 
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1 

0 


(128) 


0 0 0 

_ (72 + 7i) ^ ( 7273 ) 0 

0 ^72 (7? + 7I + 7I) ^73 (7? + ll + 7I) 0 

0 0 0 0 

It is ease to show that 7 i + 7 ^ + 73 = 0 under the condition that the grid is orthogonal at the 
boundary (Eq. (11171) ). and taking into account the Eq. fll25p . Thus, the critical matrix is: 


10 0 0 
0 C22 C23 0 

0 0 0 0 
0 0 0 0 


(129) 


with C 22 and C 23 given by Eq. (11161) . 

So, there is clearly a non-trivial incoming mode coresponding to the third right eigenvector, 
and the inflow boundary conditions are ill-posed. Looking at the critical matrix, there is a second 
incoming mode (corresponding to the fourth eigenvector), but this is actually a multiple of the first. 


because when u = +iU 


ej 


i2e2+m202 

— , tllGIl /C3 — /C4 


is ill-posed with two ill-posed modes. 


i. Hence, the initial-boundary-value problem 


6.2 Ouflow Boundary Conditions from the Second-Order Approxima¬ 
tion 

At the outflow, the generalized incoming mode may be written as 


U (^, 7 , C, t) = 


with Im{u) > 0. The wavenumber is given by: 


(130) 


k 


* 

5 


(5 + 17) (-1-S‘) 

(I5P-Z7") 


(131) 


where S* is given by Eq. (11051) . The correct square root must be taken in the definition of S* 
to ensure that if u and S* are both real then S* is positive, and if u and S* are complex then 
Im{k^) < 0. The critical matrix is actually a scalar, and has the form: 


C 


2|eP 


+ A?!702 + AiI702 + Aiie^l + A2|eCI + 2 A 1 A 2 HI 7 CI ■ 
+kiu (Aiie^l + A2|eci - - AiHIM - A2l7^) 

+0^ + AiI7l|i + 
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Under the condition that the grid is orthogonal at the boundary = |^(^| = \ri(\ = 0), (7 = 0 
leads to the same results as before: 




oj = +iU 


12Q2 + ^202 

e! 


satisfying the condition that Im{uj) > 0 , then: 


but is this case: 


r = a?02 + a202 = 




and S* = 




(132) 


(133) 




m 


r + 2i 
u 


If we substitute Eqs. fll33p into (11341) . then the result is: 


(134) 




(135) 


which is FALSE. This contradicts the assumption that there is an incoming mode satisfying the 
boundary conditions. The conclusion is that the outflow boundary condition is well-posed. 


7 Modified Boundary Conditions 

The inflow boundary conditions must be modihed to assure the well-posedness. Since the hrst three 
inflow boundary conditions require that Oi = 02 = 03 = 0 , the only nedeed condition of orthogonality 
is that between and u^. A new dehnition is proposed for V 4 : 


V4 = V4 Aimi (0 - Cx 0 0 ) -7 Asma ( 0 - 0 0 ) ( 136 ) 

The variables mi and m 2 are chosen to minimize which controls the magnitude of the 

reflection coefficient, and at the same time will produce a well posed boundary condition. The 
motivation of this approach is that the second approximation to the scalar wave equation is well- 
posed and produces fourth-order reflection. 

The new form of V 4 is: 


— ( h,l Ia,2 U,3 U,5 ) 


(137) 
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where 


^4,1 — 0 

U,2 = ix + >^iU7]^ + \2UCx - >^irrLi^y - \ 2 m 2 iz 

h ,3 = ~‘ t ~ Alt/ T]y + \2 U^y + 

= 'Cz + ^lUrjz + ^ 2 U Cz + ^27n2^x 

= l?l + + A2Fi|l 


Using Eq. (11201) . the binomial expansion of S* is: 


(138) 


= 


U 


i + U^^^j]3^|r + or) 


M 

U 


lUh^ir 

2 1 i^p 


(139) 


Thus, the product V 4 uf will be: 


V 4 uf = 


ju-m 

2\iP 


(^diA^ + y42AiA2 + 243 A 2 ) 


(140) 


where 


All = “2 ^ ^ ~ 


(141) 


^2 '^1 i^yCx ^xCy^ T ^2 {,^zVx ^x'Hz) 


(142) 


^3 — “2 “A 1^1) icT + "^2 [izCx — ixCz) (143) 

The reflection coefficient is fourth order if mi and m 2 are chosen such that V 4 U^ —)■ 0. V 4 U^ is 
a homogeneous polynomial function of Ai and A 2 ; it is zero when its coefficients are all zero. The 
problem is that there are 3 coefficients (3 equations) for 2 unknowns, mi and m 2 . 

Let’s check the consistency with the Cartesian coordinates. Using Eq. (157)) : 

Ai = -]-{u + l)-mi 

4 . = 0 

As = --{u + 1) - m 2 
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V 4 —)• 0 if ^1 = ^2 = 0 which means: 


mi = 1712 = “2 + 1) 

which is clearly consistent with the Cartesian results. The main issue here is with the term A 2 in 
equation fll40p . To close the problem, mi and m 2 need to be determined for general case. Under 
the weak assumption that ^2 ~ 0 (it must be mentioned that this assumption does not have a 
reasonable support), the values of mi and m 2 can be determined: 


mi 


m 2 


Pc^ + |g|)l')P 


i(^ + iei)icp 


The new fourth-order inflow boundary conditions will be: 


where 


■-lei 

0 

0 

0 

14 ■ 


’ ^11 

912 

913 

914 

915 

0 



0 

0 

Qt + 

921 

922 

923 

924 

925 

0 


0 

4 

0 

931 

932 

933 

934 

935 

0 


4 

4 

14 . 


. 941 

942 

943 

944 

945 . 


11 

hl 2 

hl3 

hl4 

1 - 

21 

^22 

^23 

^24 

^25 

31 

^32 

^33 

CO 

^35 

41 

^42 

^43 

/144 

-1 


9ii = 0 
gi2 = 0 
9i3 = 0 

9ia = 0 
9i5 = 0 

921 = 0 

922 = Ur]y-V^y 

923 = -Ut]^ + V^:, 

924 = 0 


(144) 


(145) 
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(146) 


and 


5'25 

^xVy ^yVx 

fi'ai 

= 0 

fl'32 

= 

fi'33 

= 0 

fl'34 

= -Uri^ + V^a: 

^735 

^xVz ^zVx 

7741 

= 0 

9A2 

= -Ur]^ + V^^ + mi^y 

7743 

= -Ur]y + V^y - 

7744 

= -Uvz + V^, 

7745 

= -c/M + fi?! 


hn 

= 

0 

hi2 

= 

0 

hi3 

= 

0 

hi4 

= 

0 

hi5 

= 

0 

^21 

= 

0 

^22 

= 

UCy - Wiy 

^23 

= 

-uCx + 

^24 

= 

0 

^25 

= 

^xCy ^yCx 

^31 

= 

0 

CO 

to 

= 


CO 

CO 

= 

0 

CO 

= 

-uCx + w^, 

^35 

= 

^xCz ^zCx 

/l41 

= 

0 

^42 

= 

~U Cx + + ^^2^2 

^43 

= 

-UCy+W^y 

/l44 

= 

-UCz + W^z- m2^x 

/l45 

= 

-I/M + 11^1^1 
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8 Summary of Dimensional Curvilinear 3D BC 

The inflow and outflow boundaries are aligned along ^ direction. 

Usefull notations: 


lel = \/(ex)' + (4)' + fe)' 

\v\ = \l iVxf + iVy? + iVzf 

ICI = \/ (Cx)^ + (4)^ + (Cz)^ 

I'Chl ^xVx T ^yVy T ^zVz (14^) 

I'CCI ~ ^xCx + ^yCy T ^zCz 

\vC\ VxCx T ^yCy T VzCz 

Tz = + 

Sa = + (6)' 


8.1 First-Order Unsteady Boundary Conditions 


The transformation to and from 1-D characteristics variables is given by the next two matrix equa¬ 
tions: 


( Cl \ 


■ -c'lei 0 0 0 lei ■ 


/ P' \ 

C 2 


0 -P^y pcix 0 0 


u' 

C3 

= 

0 -pcC 0 pcC 0 


v' 

C4 


0 pcix pciy pc^z 1^1 


w' 

^ C5 y 


0 -pc^^ -pciy -pci^ 1^1 _ 


\p' ) 




r 1 

0 

0 

1 

1 n 

1 



c2|q 

2c2g| 

25h«l 


y Cl ^ 



0 


^2 


^X 



u' 



pc^l 

2pcg|2 



C2 

v' 


0 

Cr. 

0 


^y 


C3 


pc^‘i 

2pcg|2 

2pcg|2 


w' 


0 

0 


^2 

^2 


C4 


pc^l 

m\i? 

2pcg|2 


\v' ) 


0 

0 

0 

W\ 

M - 


\ C5 / 


(149) 


( 150 ) 


where ci C 2 C 3 C 4 and C 5 are the amplitudes of the five characteristics waves. 
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At the inflow boundary, the correct unsteady, non-reflecting boundary conditions for a subsonic 
flow are: 


while at the outflow boundary: 


/ Cl \ 


C2 

C3 


= 0 



(151) 


C5 — 0 

8.2 Approximate, Quasi-3D, Unsteady Boundary Conditions 

At the inflow, the boundary conditions in terms of flow perturbations are: 


where 


c^lel 

0 

0 

0 

lel ■ 


0 

-P^y 

pcix 

0 

0 

Qt + 

0 

-pC^z 

0 

pcix 

0 

0 

pcix 

P^y 

pciz 

lel . 

r hii 


h-21 

hsi 

hii 


9ii 

fi'12 

fl'is 

914 

921 

922 

923 

924 

931 

932 

933 

934 

. fi'41 

942 

943 

944 

hi2 

hl3 

hi4 

hl5 

^22 

^23 

^24 

^25 

^32 

^33 

^34 

^35 

^42 

^43 

^44 

^45 


Q. 


Qc = o 


Q = ( p' u' v' w' p' )" 


and at the outflow, they are: 


[ 0 -pc^x -pciy -pci^ 1^1 ] Qt + [ 5'51 5'52 5'53 5'54 5'55 ] Qr, 

+ [ hsi /l52 /i53 ^54 hss ] = 0 


where 


(152) 


(153) 


(154) 


(155) 
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and 


gii = 0 

912 = 0 

913 = 0 

914 = 0 

915 = 0 

921 = 0 

922 = Uriy-V^y 

923 = -Urjx + Vi^ 

924 = 0 

925 ^xVy ^yVx 

931 = 0 

932 = (156) 

933 = 0 

934 = -Urjx + Vi^ 

935 ^x^z ^z9x 

941 = 0 

942 = -Urj^ + V^x 

943 U9y~^ ^ 

944 = -Urjz + Vi^ 

945 — -U-^ + V\Q 

951 = 0 

952 = 

953 = Urfy-Viy 

954 = Ur]^-V^z 

te = -C7J||1 + F|^| 


hii = 0 

hi2 = 0 

/ll3 = 0 

hi4 = 0 

hi5 = 0 
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h2i — 0 

h22 = UCy-W^ 
h23 = -UC, + W^, 

/i24 = 0 

^25 = ^xCy ^yCx 

hzi = 0 

h32 = UC.-Wi, 

^33 = 0 

h3A = -UCx + Wi:, 

^35 ^xCz ^zCx 

h^i = 0 

h42 = -UCx+W^, 

h^3 = -UCy + my 

h^A = —UC,z + W^z 

hit. = -C7^+W|«| 

^51 = 0 

h^2 = UCx-Wi^ 
h3 = UCy-my 
hA = UC.-Wiz 

htt = -C7^+W|«| 

At the inflow, the boundary conditions in terms of characteristic variables are: 


/ Cl \ 


C 2 

+ 

C 3 


\ C 4 y 

hii 


^21 

+ 

CO 




and at the outflow, they are: 




f Cl \ 

9 u 912 9 i 3 9 ii fi'is 

fl '21 922 923 92 A 925 

d 

C 2 

931 932 933 934 935 

dr] 

C 3 

941 942 943 944 945 . 


C 4 

\C 5 ) 


hi 2 

hl 3 

hiA 

hl 5 


( Cl \ 

C 2 

^22 

^23 

^24 

^25 

d 

^32 

^33 

^34 

CO 


C 3 

^42 

^43 

^44 

^45 J 


C 4 

V C 5 / 


(157) 


(158) 
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where 



^ Cl ^ 


C2 

d 


dr] 

C 3 


C 4 


\ C 5 ) 

( Cl 

\ 


d 

hsi /l52 /i53 /l54 /l55 — 

L J (9C 


C2 
C3 
C4 

VC5/ 


= 0 


911 

912 

913 

914 

915 

921 

922 

923 

924 

925 

931 

932 

933 

934 


0 

0 

0 

0 

0 

0 


^2 + ^xVx) + V 

{UVy-ny) 

0 


(Uv, - r&) 

U , , — 

“ 1 “ ^xVx) + V 


V 2|ep J 


i^xVz 


^zVx) 


(159) 


(160) 
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and 


fl'35 

fi'41 

9i2 


953 

fl'54 

955 


) i^xVz - ^zVx) 


= 0 


^2 (^y9x ^xVy) 

94:3 ~ “172 i^zVx CxVz) 




944 — 

945 = 

951 = 

952 = 


_t/IM + V 
lep 


0 

0 


[/ 


942 ^2 (^y9x CxVy) 

fi'43 = i^zVx ^xVz) 


= 0 


jj\^v\ .TT 
944 - ^ 


hii 

hi2 

hi3 

hiA 

hi5 

^21 

^22 

^23 

^24 

^25 

^31 

^32 
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0 

0 

0 

0 

0 


^2 (^y^y ^xCx) + W 

m + u 

2|{P_ 

0 


i^xCy ^yCx) 
i^xCy ^yCx) 


-^4 {UC. - wQ 


( 161 ) 
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^33 

hzA 

^35 

hil 

hA2 

hiS 

hiA 

/l45 

^51 

^52 

^53 

/l54 

^55 


+ Ca:Cx) + 


c|el + t/ 
2|eP_ 
c\i\-U 
2|eP 


i^^xCz ^zCx) 
i^xCz ^zCx) 


i^yCx ^xCy 


i^zCx ^xCz) 


0 

u 

n 

-ztM+tt 

0 

0 

—/l 42 — i^yCx ~ ^xCy) 

^43 (CzCa; ^xCz) 

^3 

0 

/i44 = -r/^+W 


lel^ 


8.3 Modified Boundary Conditions 

Under the assumption that 2 I 2 ~ 0 (?; this assumption need to be revised), the values of mi and m 2 
could be determined: 


mi 


m 2 


Hc'+kl)l>)|- 

^yVx ^xVy 

H^+iei)lcP 

^zCx ^xCz 


(162) 


The new fourth-order (?) inflow boundary conditions in terms of flow perturbations will be: 


■ -c^lel 

0 

0 

0 

lel ■ 


' 911 

912 

913 

9ia 

915 

0 

-pc^y 

pcix 

0 

0 

Qi + 

921 

922 

923 

92A 

925 

0 

-pC^z 

0 

pcix 

0 

931 

932 

933 

93A 

935 

0 

pcix 

P^y 

pciz 

lel . 


. 9ai 

9A2 

9A3 

9aa 

9A5 . 
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Qc = o 


(163) 


where 


and 


■ h 

h 

h 

h 


11 

hl2 

CO 

hi4 

hl5 

21 

^22 

^23 

^24 

h25 

31 

^32 

CO 

CO 

^34 

^35 

41 

^42 

CO 

/l44 

/l45 _ 


gn = 0 
9i2 = 0 
9n = 0 
gu = 0 
fi'is = 0 
921 = 0 

922 = Urjy - V^y 

923 = —Urjx + V^x 

924 = 0 

925 — ^xVy ^yVx 

931 = 0 

932 = Ut]^ - ( 164 ) 

933 = 0 

5-34 = -Urix + V^x 

935 ^xVz ^zVx 

941 = 0 

942 = -Ur]x + V^x + mi^y 

943 = -Urjy + V^y-mi^x 

944 = -Urj^ + V^z 

945 = 


hii = 0 

hi2 = 0 
hi3 = 0 

hi4 = 0 

hi5 = 0 
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h2i — 0 

h22 = UCy-W^ 
h23 = -UC, + W^, 

/i24 = 0 

^25 = ^xCy ^yCx 

hsi = 0 

h32 = UC.-Wi, (165) 

^33 = 0 

/^34 = -UCx + Wi, 

^35 ^xCz ^zCx 

/i41 = 0 

^42 = ~U(^x + W + TT^2 ^z 

h43 = -UCy + my 

/i44 = -UCz + W^^- m2ix 

fc« = -c/^ + w"l«l 

At the inflow, the new boundary conditions in terms of characteristic variables are: 


Wt 


f Cl \ 


C2 

+ 

C3 


\ C4 y 

hii 


^21 

+ 

^31 


hil 


9ii 

921 

931 

9ai 

hi2 

^22 

^32 

^42 




f Cl \ 

9l2 9l3 9lA 915 

922 923 924 925 

d 

C2 

932 933 934 935 

dr] 

C3 

942 943 944 945 . 


C 4 

V C5 y 


hi3 

hi4 

hl5 


y Cl ^ 

C 2 

to 

CO 

^24 

^25 

d 

^33 

^34 

^35 


C3 

CO 

^44 

/l45 


C 4 

\C5 J 


and at the outflow, they are: 


(166) 
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/ Cl \ 


where 


+ \_ fl'51 9^2 953 954: 955 ] ^ 


+ /isi h. 


52 ^53 /l54 /l55 1 — 

J <9C 


C2 
C3 

C4 

\C5/ 

/ Cl \ 

C2 
C3 
C4 

\C5J 


= 0 


911 

912 

913 

914 

915 
921 

922 

923 

924 

925 

931 

932 

933 

934 
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0 

0 
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0 

0 


^2 + V 

0 


(U’h - r&) 

U , , — 

+ ix9x) + V 


V 2|ep J 


i^xVz 


^zVx) 


(167) 


(168) 
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- izTlx) 


i^yVx - ixVy) - mi 

U , , U 

{^zVx ~ ^xVz) ~ mi-^^y^z 

_^IM + 7 
lep 


^2 i^yVx CxVy) 


i^zVx ^xVz 


.^IM + y 

ler 


^2 (^yQ ^xCx) + W 


(u, - e,cj 

(«.c, - ?»c.) 


m. - »'(.) 


(169) 
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^33 

^34 

^35 

/l41 

^42 

^43 

/l44 

/l45 

^52 

^53 

/l54 

^55 


■4 (t-t+ &« + »' 

(&C. - U.) 


2|«l _ 
SKI - U 
2|?P 


i^xCz ^zCx) 


0 




^2 ^xCy) ^2 ,-r,2 'CyG 


('CzCx ^xCz) ^2 


-uM + 

0 

0 

('Cj/Ca; ~ ^xCy) 
(CzCx CxCz) 

^3 

0 

-uMl 

lep 
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